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Abstract 

Parameter determination during an injectivity test relies on analytical models that fully foresee how pressure 
behaves when water is injected in the reservoir. Existing analytical models are able to depict both injection 
and falloff periods in single-layer reservoirs, but lack to describe the falloff period in multilayer systems. 
Therefore, this work attempts to extend the existing single-layer falloff formulation to multilayer reservoirs. 
The accuracy of the proposed solution was assessed by comparison with a numerical simulator. The suggested 
model was also used to determine the reservoir equivalent permeability. 
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1. Introduction 

The injectivity test is a procedure used to collect 
information over a petroleum reservoir by injecting 
water into the reservoir. Several reservoir features, 
such as equivalent permeability and outer boundary 
condition, might be inferred from the pressure 
response measured during the test. 

Injectivity test consists of two different stages: 
the injection period and the falloff period. During 
the former, occurs the water injection into the rock 
formation. The latter stage is marked by the well 
shut-in and, hence, a zero-flow pulse propagates 
along the reservoir. 

Over the past years, accomplishments have 
been made regarding the pressure behavior in 
multilayer reservoirs under single-phase flow and 
injectivity tests in single-layer reservoirs. However, 
an analytical solution for pressure behavior in 
multilayer reservoirs is well known just during the 
injection period. 

Therefore, this work attempts to develop a 
new analytical model for the falloff period in 
multilayer reservoirs. The suggested formulation 
was reached by combining the solution for injection 
period in multilayer reservoirs with the existing 
model for falloff period in single-layer reservoirs. 
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The developed formulation was, then, applied to 
estimate the reservoir equivalent permeability from 
pressure data. 

2. Previous Achievements 

Historically, pressure behavior in petroleum 
reservoirs has been modeled under two distinct 
lines of work. On one hand, analytical models 
for pressure response in multilayer reservoirs under 
single-phase flow have been developed. In parallel, 
the study on injectivity tests in single-layer 
reservoirs has been conducted. 

2.1. Modeling Single-Phase Flow in Multilayer 
Reservoirs 

Single-phase flow in multilayer reservoirs has 
been studied in the context of production tests. 
However, the same equations also describe a 
theoretical injectivity test where the reservoir 
existing oil and the injected fluid present the same 
properties. 

Lefkovits et al. [1] provided an analytical model 
for production tests in commingled systems. Their 
model, which has been served as basis for many 
studies further developed, considers that layer 
properties may be different in each layer. 

Cobb et al. [2] conducted a comparison 
between three techniques for buildup data 
analysis in multilayer commingled reservoirs 
with equal layer thickness: Muskat, Horner and 
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Miller-Dyes-Hutchinson methods. They have 
shown that all three mentioned methods might be 
used to estimate the reservoir flow capacity. 

Single-phase flow in multilayer reservoirs with 
two layers of distinct thickness was studied by 
Raghavan et al. [3]. A means for obtaining 
individual layer permeability was proposed, 
assuming that thickness ratio is known. 

Gao [4] developed approximated solutions for 
pressure and flow-rate in each layer. The suggested 
formulation is able to provide estimates for flow 
capacity in each layer. 

A procedure for individual layer properties 
evaluation was reached by Ehlig-Economides and 
Joseph [5]. The proposed calculation relies on the 
analysis of pressure and layer flow-rate data. 

2.2. Modeling Injectivity Tests in Single-Layer 
Reservoirs 

The study on injectivity tests started to drain 
greater attention in the 80’s. Abbaszadeh and 
Kamal [6] proposed the waterfront saturation 
profile depicted by Buckley-Leverett theory [7] 
may be discretized into a series of banks. Thus, 
pressure behavior during an injectivity test would 
be equivalent to the one of a composite reservoir. 
They also proved that, under certain conditions, 
waterfront remains stationary during falloff. 

Bratvold and Horne [8] evaluated the influence 
of temperature on fluid mobility and on saturation 
gradients. Their formulation was also based 
on discretizing the saturation profile foreseen by 
Buckley-Lever ett theory to compute the pressure 
response. 

A theory for both single-phase and multiphase 
flow in radially heterogeneous reservoirs was 
developed by Thompson and Reynolds [9]. They 
showed that pressure derivative may be understood 
as a weighted average of the permeabilities along 
the reservoir. The weighting factor is a function of 
flow-rate and mobility gradients. 

In a subsequent work, Banerjee et al. [10] 
modeled pressure behavior during injectivity tests 
in radially heterogeneous reservoirs. The same 
suggested formulation also provides a means for 
computing the mechanical formation damage in 
homogeneous reservoirs. 

A general theory for the injection period was 
presented by Peres et al. [11]. They showed that 
their model reduces to the formulation proposed 
by previous authors, provided that the correct 


assumptions are made. Their solution is also 
capable of describing restricted entry wells. 

The falloff solution was reached by Peres et 
al. [12], following their own work [11]. The 
superposition principle was used to describe the 
flow-rate history during falloff as the sum of two 
flow-rates with different sign being applied at the 
same point of the reservoir. 

2.3. Modeling Injectivity Tests In Multilayer 
Reservoirs 

The first work regarding injectivity tests in 
multilayer reservoirs was presented by Barreto et 
al. [13]. They applied Darcy’s law to one 
given layer, so that an expression for the pressure 
change during the injection period in this layer is 
obtained. Thereby, layer flow-rate may be written 
as a function of that individual layer pressure 
change. Wellbore pressure during injection is, then, 
computed by summing up all layer flow-rates. Their 
solution assumes no crossflow and radially infinite 
reservoir. 

As disclosed in this overview, existing analytical 
models lack to describe falloff period in multilayer 
systems. Using the formulation reached by Barreto 
et al. [13] as basis, this work will attempt to extend 
the analytical model for falloff period in single-layer 
reservoirs [12] to multilayer systems as well. To the 
author’s best knowledge, the developed formulation 
has not yet been presented before. 

3. Reservoir Model 

The analytical model developed was based on 
two fundamental principles: the analytical model 
for falloff period in single-layer reservoirs [12], and 
the existing formulation during the injection period 
in multilayer systems [13]. All calculations assume 
that a consistent set of units is used and that 
the reservoir is subject to the following simplifying 
hypothesis: 

• At the instant t = 0, reservoir is in equilibrium; 
z.e., pressure is the same in all layers; 

• Homogeneous and isotropic reservoir, with 
infinite extension; 

• Water and oil are assumed to be imiscible, 
slightly compressible fluids with constant 
viscosity /q 

• Flow is isothermal; 
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• Rock formation presents a low and constant 
compressibility; 

• Wellbore fully penetrates all layers and injects 
at constant flow-rate qwy, 

• Gravitational and capillary forces are 
neglected; 

• There is no wellbore storage; 

• It is considered a commingled system. 

Figure 1 displays the considered reservoir model. 
Reservoir may be composed by an arbitrary number 
of layers, which may present distinct formation 
damage zones, or even no formation damage at 
all. Oil properties, relative permeability curves and 
rock porosity were assumed to be the same in all 
layers during the calculations. 
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Figure 1: Reservoir Model. 


4. Modeling Two-Phase in Multilayer 
Reservoirs 

Darcy’s law, which models flow through porous 
media, is the starting point to compute the pressure 
variation in the reservoir during an injectivity test 
[ 10 ]. 

Since the reservoir is composed of multiple layers, 
Darcy’s law is applied in one given layer j [13]. So, 
considering the reservoir geometry and the model 
hypothesis, oil and water flow-rates are given by: 


/ . 2irhjknk ro dPj 

Qoj (r,t) = — - 


Mo 


dr 


and: 


(i) 


7 x 27r hnknk rw dPj 

Qwj ( r ,t) = — - 3 3 rw 3 


k J W 


dr 


( 2 ) 


Equations (1) and (2) are reached by assuming 
that relative permeability curves and oil viscosity 
are the same in all layers. This assumption is made 
only for a matter of convenience. Both equations 
can be easily modified if fluid properties are distinct 
for each layer. 

Mass conservation inside the reservoir states 
that total flow-rate at any layer of a commingled 
reservoir is computed as: 


qj (r, t) = q wj (r, t) + q Q j (r, t) = 


0 7 7 / ^ro . k T W 

— ZTTfljkj -t- 

V Mo Mw 


dP±_ 

dr 


( 3 ) 


At this point, it is convenient to write equation 
(3) in terms of fluid mobility. This property 
measures the easiness of one phase to move 
throughout the porous media. It is defined as the 
ratio between the phase’s relative permeability and 
its viscosity [9]: 


A f (S w ) = (4 ) 

Mi- 

Total mobility of the system in study is given by 

[ 12 ]: 


Xt(S w ) = X 0 (S W ) + X W (S W ) (5) 

During an injectivity test, a region with high 
water saturation is formed around the wellbore. As 
injection goes on, this flooded area grows. Hence, 
water and oil saturations are functions of time and 
radius. Then, relative permeabilities and mobilities 
will also depend on time and radius. 

Thereby, Darcy’s law may be rearranged, so that 
pressure change at any point of the reservoir is 
expressed as a function of the radius: 


dR, 

dr 


Qj(r,t) f _1_\ 1 

2 nkjhj \X 0 (r,t) + X w (r,t) J r 
= Qj(r,t) 1 1 

2tt kjhj At(V, t) r 


(6) 
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Thus, integrating both sides over the radius, from 
the wellbore onwards [9]: 


A Pwfj (t) 


oo 

1 f dr 

2tt kjhj J \r{r,t) r 

r w 


( 7 ) 


where the wellbore radius is denoted by r w . 

Beyond the waterfront, only oil flows. On the 
other hand, behind the waterfront, a two-phase 
flow occurs [6, 8]. Hence, it is logical to split 
the right side of equation (7) in two integrals, one 
whose integration limits go from the wellbore to 
the waterfront and another whose limits go from 
the waterfront onwards: 


A Pwfj (t) 
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In the second integral of equation (8), total 
flow-rate was replaced by oil flow-rate at irreducible 
water saturation (g OJ ), and total mobility was 
replaced by the oil mobility at irreducible water 
saturation (A 0 j). This comes from the fact that 
oil is the only fluid flowing beyond the waterfront 
[6, 8]. Furthermore, oil and water remain at 
their respective connate saturations beyond the 
waterfront, since injected water has not yet reached 
this region. Thus, oil properties must be evaluated 
at its endpoint saturation. Waterfront radius, 
denoted by rpj (t), is estimated by Buckley-Leverett 
theory [7]. 


5. Proposed Formulation 


Equation (8) concerns a two-phase flow in a 
reservoir subject to the assumptions made in 
section 3. It applies for both injetion and falloff 
periods. The solution for the injection period was 
derived from equation (8) by Barreto et. al [13]. 

To reach the formulation for the falloff period, 
additional assumptions must be made, regarding 
the physical phenomena that occur after the 
wellbore is shut. 


5.1. Pressure Behavior During Falloff 

It is known that the mobility profile remains 
practically unchanged during falloff [6, 8]. 

This hypothesis plays an essential role to the 
computation of wellbore pressure during falloff. If 
the waterfront were transient, flow-rate throughout 
the reservoir would also be transient . Pressure 
behavior, then, could only be solved numerically. 

Thus, in equation (8), total mobility profile and 
waterfront radius depend only on the radius and on 
the injection time t p : 


r Fj ( t p ) 
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+ 


OO 

— f 

2irkjhj J 


r F, ( t p ) 


Qsjjr, At) dr 
A r(r,t p ) r 

dosAr, A t) dr 
A „ r 


( 9 ) 


where the subscript s denotes the falloff 
period. Adding and subtracting the term 
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equation (9): 
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in the right side of 


AP WSj (At) = 
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Equation (10) shows that the pressure variation 
at the wellbore may be understood as the sum of 
two terms: one owing to the mobility differences 
between oil and water (AP\ Sj ), and another related 
to the single-phase oil flow (A P OSj ): 


AP WSj {At) = AP Xsj (At) + AP OSj {At) (11) 

where: 
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and: 


AP OS j(At) = 


OO 

— / 

2nkjhj J 


qosjjr. At) dr 

A„ r 


(13) 


By hypothesis, there is no formation crossflow 
and pressure is the same in all layers, apart from 
the hydrostatic effect. This means that, after the 
well is shut, the system is in equilibrium and there is 
no flow between layers, neither through the vertical 
boundaries between adjacent layers nor through the 
wellbore. This assumption is not only crucial to 
the development of the falloff formulation, but also 
significantly strong. 

Assuming that there is no flow between layers 
means, in other words, to state that vertical 
flow between layers is so small compared to the 
horizontal flow along the reservoir, that it may be 
neglected. This might not be true during falloff if 
layer properties are remarkably different [14]. 

From the negligible vertical flow hypothesis, 
follows that the zero-rate pulse equally propagates 
along all layers during falloff. Then, total and oil 
flow-rate may be expressed as functions of layer 
flow-rate fraction, which is the ratio between the 
flow-rate in a given layer j evaluated at one given 
radius and one given instant of time, and the 
injected flow-rate: 


9D ,M) = ^ (14) 

Qinj 

To reach the expression for pressure change 
during falloff, layer flow-rate fraction within the 
flooded region is assumed to remain constant after 
the well is shut, and equal to the flow-rate fraction 
at t = t p . Then, total and oil flow-rates in layer j 
may be expressed as: 


q sj (r,At) = q Dp -q s (r,At); 

and: 

q OS j(r,At) = q Dpj q os (r,At), 


(15) 

(16) 


where qu Pj is the layer flow-rate fraction just before 
shut-in, defined as: 


QD Pj = 


qj(t — t p ) 
Qinj 


(17) 


AP WSj (At) = AP OSj (At)+ 

r Fj (tp) 

VDpj [ (0s_ _ qos\ dr 

2 nkjhj J \X T A Jr K > 

r w 

It is known that the term A P os j(t) is, in fact, 
the same for all layers [4, 5, 13]. From now on, this 
term will be referred to as A P os (t). Furthermore, 
the considered reservoir model states that pressure 
change is the same for all layers. So, the pressure 
measured at the well bottom is equal to the pressure 
at any individual layer, apart from the hydrostatic 
effect: 


AP WSl (At) = ... = AP WSn (At) = AP ws {At) (19) 

Besides, the term multiplying the layer flow-rate 
fraction will be denoted as a weight variable 


R 3 (At) 


1 

2irkjhj 


r F j (t P ) 

1 
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A T 


qos \ jr 
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( 20 ) 


Thus, an expression to compute qD Pj is obtained 
from equation (18): 


QD Pj 


AP ws (At) - AP os (At) 
Rj(At) 


( 21 ) 


At any time step, in particular for t = t p , the 
definition of layer flow-rate fraction ensures that: 


ft ^ 76 11 

T dDj = - qj =^-l =>■ T dDpj = 1 (22) 

i=l qinj 3=1 j =1 

Hence, from equations (21) and (22): 


■ (AP ws {At) - AP os (At)\ 

“1V RA At) ) 

= (AP ws (At) - AP os (At)) (]^y) (23) 


Thus, equation (10) may be rewritten as: 


Rearranging equation (23): 
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A P ws (At) = AP os (At) 


y - 


(24) 


As described in equation (24), pressure change 
in multilayer reservoirs during falloff may be 
understood as the sum of the single-phase 
contribution with one term that encompasses the 
mobility differences between water and oil, in an 
analogous way to the injection period. 

In single-layer reservoirs, equation (24) reduces to 
the analytical model for falloff period in single-layer 
reservoirs proposed by Peres et al. [12]. This 
ensures that the suggested formulation is able 
to describe pressure behavior during falloff in 
reservoirs with an arbitrary number of layers. 


5.2. Flow-Rate Approximation 

The comprehension of the procedure described 
in section 5.1 shows that an essential step consists 
of determining flow-rate at any point along the 
reservoir. Thereby, equation (24) may be applied 
to compute the pressure change during falloff. 

One central question is how to estimate 
total flow-rate behind the waterfront. Two 
straightforward and rough approximations come 
from considering that total flow-rate equals either 
to oil or to water flow-rate, calculated at their 
respective endpoint saturations: 


q s (r, At) = q os {r, At) 

(25) 

q s (r,At) = q ws (r, At) 

(26) 


A more accurate approximation is obtained by 
taking an average between these two flow-rates, 
weighted by water and oil endpoint saturations: 


An estimation for water and oil flow-rates may be 
derived from the dimensionless version of Darcy’s 
law: 


q D (r,At) = -r D dPr £ r,At \ (29) 
or D 

where gfo, td and Pd are the dimensionless 
flow-rate, radius, and pressure change, defined as: 


q 

qD = —; 

Qinj 

(30) 

r 

r D = —; 

(31) 

khX f AP 

Qinj Ff 

(32) 


The subscript / indicates the flowing phase. 
Therefore, flow-rate may be easily computed once 
the pressure change is well defined. For radially 
infinite flow, a simple and quite accurate estimate 
for the flow-rate is provided by the source line 
solution [12]: 


qD=exp {y) <33 > 

where the dimensionless time to is computed as 
follows: 


I'D 


A fkt 

<pc t rl 


(34) 


At this point, another fundamental assumption 
is made. In multilayer reservoirs, it was assumed 
that, if the dimensionless parameters are calculated 
using the reservoir’s total thickness and equivalent 
permeability, then the source line approximation 
yields the sum of oil (or water) flow-rate in all 
layers. Total thickness and equivalent permeability 
are defined as [1, 2]: 


</.-('■• At) = 


Kqos(r, At) A wQwsfri At) 

A 0 “t“ \yj 


(27) 


Finally, approximation (27) is improved by 
accounting for the mobility profile along the flooded 
area: 


q s (r, At) 


yQqosip At) + A w (r)q ws (r, At) 
A c (r) + X w (r) 


n i n 

hr = hj ; k eq = —— kjhj (35) 

3=1 T 3=1 

Definitions for dimensionless radius and flow-rate 
remain the same as in the single-layer case. 

During falloff, pressure change is computed by 
applying the superposition principle, as depicted by 
Peres et al. [12]: 


(28) 
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Oil flow-rate at a given radius is evaluated 
by applying equation (36). Oil properties 
should be used to compute the dimensionless 
parameters. On the other hand, water flow-rate 
is calculated by using water properties to compute 
the dimensionless variables. 

5.3. Adjustments Due to the Skin Effect 

So far, all calculations assume the reservoir does 
not present any formation damage. Nevertheless, 
perforation and completion of the wellbore may 
create a region with modified permeability ( k s ki n ) 
around it. This phenomena, known as skin effect, 
results in an additional pressure rise, or drop, 
depending on the value k s ki n - 

With the purpose of accounting for such effect, 
the Rj coefficient used in the calculation of AP\ 
must be adjusted whenever the reservoir presents 
formation damage. The correction depends on the 
waterfront position, compared to the skin radius 
(rskin)' While the flooded region is within the 
damaged zone: 


r Fj (tp) 



7 *, 


{2,7) 

Otherwise, if the waterfront has overcome the 
damaged zone: 


rFAt p ) 



6. Results and Discussion 

To evaluate the accuracy of the solution proposed 
in section 5, a set of cases was run on a commercial 
based flow simulator. Oil model used was blackoil. 

Pressure variation is expected to change in time 
as two semilog straight lines, one associated to the 
oil properties and another governed by the water 
properties, separated by a transition period [9, 10]. 

Therefore, in order to properly interpret the 
results of an injectivity test, one must analyze 


krel 



Figure 2: Relative Permeability Curves 


not only the pressure data, but also the pressure 
derivative with respect to the logarithm of time 
(from now on, referred to as pressure derivative). 
Each semilog straight line is identified by a period 
when pressure derivative is constant, and reflects 
the properties of one distinct fluid. 

6.1. Input Parameters 

For all cases, it was considered a 4 days (96 hours) 
injection period followed by a 4 days falloff period, 
as in a typical injectivity test. The 4 days injection 
time is short enough so that infinite radial regime is 
still acting [15]. Injection flow-rate, wellbore radius 
and other rock and fluid properties may be seen in 
table 1. 

Oil and water relative permeability curves, which 
were assumed to be the same in every layer, 
are displayed in figure 2. Water properties 
were evaluated at a typical reservoir temperature 
(between 50 and 60 °C). 

For all cases, total reservoir thickness was set as 
30 m, to ensure that the hypothesis of negligible 
gravitational effects is valid. 

Depending on the relative easiness of water 
and oil to move through the reservoir, the 
two-phase flow may be labeled as either favorable 
or unfavorable to water displacement. One useful 
parameter to understand which fluid will moves 
more easily is the endpoint mobility ratio M: 

M=A (39) 

Flows favorable to water displacement imply that 
M > 1, while mobility ratios lower than 1 denote 
the flow is unfavorable to water displacement. 
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Qinj (m 3 /d) 

r w (m) 

4> 

c r (kgf/cm 2 ) 1 

l^w (cP) 

c w (kgf/cm 2 ) 1 

c 0 (kgf/cm 2 ) 1 

500 

0.108 

0.32 

8.0xl0 -5 

0.52 

i.i4xi(r 4 

4.04xlCT 5 


Table 1: Rock and Fluid Properties 


For each reservoir configuration, a pair of values 
for oil viscosity was chosen such that one would 
result in a flow favorable to water (/i 0 = 5.10 cP) 
and the other would present a mobility ratio lower 
than one (/i G = 1.00 cP). Table 2 shows the reservoir 
settings for each case. 

6.2. Comparison Between Analytical Model and 
Numerical Simulation 

Results for cases depicted in table 2 are displayed 
in figures 3 to 8. Analytical solution and numerical 
data showed a good agreement for all cases. 

Cases A1 and A2 show the pressure response 
in a reservoir without formation damage and with 
different permeabilities in each layer. Thickness is 
the same in all layers. Difference between individual 
layer permeabilities is intended to assess whether or 
not the hypothesis of no flow between layers during 
falloff is valid. 

The pressure behavior in a two-layer reservoir 
with two distinct damaged regions is displayed in 
cases B1 and B2. This case aims to evaluate 
the influence of distinct skin factors in the overall 
pressure behavior. 

Cases Cl and C2 present five distinct layers. 
Damaged zone properties imply that skin factor 
is the same in every layer. Again, the hypothesis 
of negligible flow between layers during falloff is 
tested. 



6.3. Results Analysis and Discussion 

It is noticeable that pressure derivative starts at 
a level during early time and, after a transition 
period, stabilizes at a different level at long time. 
This is a typical feature of a two-phase flow. 

At early falloff times, pressure is governed by 
water properties. As the zero-rate pulse reaches the 
higher oil saturations regions, pressure derivative 
starts to reflect the properties of this fluid [16]. 
For cases whose flow is favorable water, the 
highest pressure derivative level is associated to 
oil. Similarly, for cases with M < 1, the highest 
derivative level is related to water. 

Differences between simulated and analytical 
data are related to the fact that numerical 
simulation faces convergence issues due to the 
existence of shocks in the solution [13]. Thus, the 
transition region between water and oil properties 
level may be delayed in the numerical simulation. 

Early falloff times present some divergences 
between numerical and analytical data, especially 
in cases B1 and B2. The main reason for that relies 
on usage of the source line solution to determine 
total and oil flow-rates. Such approximation is less 
accurate in the region closest to the wellbore. As 
the zero-rate pulse develops, this error becomes less 
significant. 

Moreover, the choice of the total flow-rate 
approximation is another cause of deviations 
between numerical and analytical data during 
early falloff times. Figure 9 compares the 



Figure 3: Case A1 


Figure 4: Case A2 
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Case 

Layers 

k (mD) 

h (m) 

kskin (mD) 

r s kir 

i M 

M 



1500 






Al 

3 

500 

10 (all) 

- 


- 

3.09 



1000 








1500 






A2 

3 

500 

10 (all) 

- 


- 

0.61 



1000 






B1 

2 

1000 (all) 

15 (all) 

500 

100 

0.5 

(all) 

3.09 

B2 

2 

1000 (all) 

15 (all) 

500 

100 

0.5 

(all) 

0.61 



1000 


500 






800 


400 




Cl 

5 

1200 

6 (all) 

600 

0.5 

(all) 

3.09 



900 


450 






1400 


700 






1000 


500 






800 


400 




C2 

5 

1200 

6 (all) 

600 

0.5 

(all) 

0.61 



900 


450 






1400 


700 





Table 2: Tested Cases 
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wf Analytical 
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^ Oil Properties Level 
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At (h) 


Figure 5: Case B1 


Figure 6: Case B2 


results obtained by the analytical model using 
approximations 1 to 4 (equations (25) to (28)). 
Reservoir and fluids settings are the same as case 
Al. 

When approximation 1 (equation (25)) is used, 
the overall derivative behavior is consistent with 
a two-phase flow. However, approximations 2, 3 
and 4 (equations (26) to (28)) yielded a pressure 
derivative profile incompatible with the physical 
reality, presenting an oscillatory region before the 
oil properties level is reached. 

Such unexpected fact is possibly related to the 


development of the formulation for the injection 
period. The model proposed by Barreto et al. 
[13] is achieved under the hypothesis that the 
flooded zone is within the steady state region. 
Hence, inside the waterfront, total flow-rate and 
oil flow-rate at irreducible water saturation are 
numerically equal during the injection period. 
Perhaps the results displayed in figure 9 reflect this 
assumption, suggesting that, during falloff, total 
flow-rate should be estimated by the oil flow-rate. 

Assuming the reservoir is in hydrostatic 
equilibrium after shut-in is a quite strong 
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Figure 7: Case Cl 


Figure 8: Case C2 


hypothesis, specially when layer permeabilities and 
damaged zones are very different. For these cases, 
flow between layers through the wellbore may occur 
right after the well is shut, due to the differences 
between layer flow capacities [4, 14]. Neglecting this 
effect is a possible cause of the divergences between 
analytical and numerical pressure derivatives at 
early falloff time (At < 0.001 h). 

Besides the flow-rate approximation, another 
possible explanation for this lower derivative level 
is the accuracy of the Bessel functions used to 
evaluate A P os term. The computation of A P os 
involves modified Bessel functions Kq and K\ [13]. 
For very small time arguments, the imprecision 
associated to the computational calculation of those 
functions becomes more relevant. 

Thus, initial falloff times are subject to an error 
that is inherent to the flow-rate approximation 
and computational issues, but not related to the 
hypothesis required to reach equation (24). Despite 
that, as the zero-rate pulse propagates through 
the reservoir, analytical solution and numerical 
data rapidly converge, as the zero-pulse propagates 
throughout the reservoir. 

Divergences between analytical solution and 
numerical simulated data are also related to the 
numerical calculation procedure. Flow simulator 
assumes that water is injected by a zero-radius 
wellbore, that is, well is considered to be a source 
line. Moreover, the pressure is evaluated as the 
average pressure in each grid block. This implies 
the evaluated pressure is influenced by the radial 
steps in the simulation grid [17]. Although such 
factors are inherent to the numerical simulation, 
they are not an issue for the analytical model. 

Yet, the errors associated to the flow-rates 


computation are more relevant in cases with 
formation damage. The higher flow-rate values 
estimated by the source line approximation imply 
in slightly higher values of A P\ when total flow-rate 
is estimated through equation (25). In cases with 
M > 1, the pressure term associated to the 

mobility differences is negative. Hence, during 
initial falloff times, this reflects in lower pressure 
derivative values. In cases whose flow is unfavorable 
to water, this error is hidden, since the A P\ 
term is already positive. As the zero-rate pulse 
propagates throughout the reservoir, this effect 
gradually vanishes and analytical and simulated 
data converge. 

6-4- Determining the Reservoir Equivalent 
Permeability 

As mentioned in section 6.3, two flatten regions 
are seen in the pressure derivative profile: one 
defined by the oil mobility and another related 



0,001 0,01 0,1 1 10 100 
At(h) 

Figure 9: Falloff Data for Case A1 Using Different Flow-Rate 
Approximations 
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Case 

k eq (mD) 

Phase 


k eq Calc. (mD) 

Error (%) 

Al 

1000 

Oil 

1.483 

1010 

1.0 

Water 

0.508 

955 

-4.5 

A2 

1000 

Oil 

0.295 

995 

-0.5 

Water 

0.522 

929 

-7.1 

Bl 

1000 

Oil 

1.486 

1008 

0.8 

Water 

0.351 

1382 

38.2 

B2 

1000 

Oil 

0.295 

995 

-0.5 

Water 

0.536 

905 

-9.5 

Cl 

1060 

Oil 

1.399 

1070 

1.0 

Water 

0.505 

960 

-9.4 

C2 

1060 

Oil 

0.278 

1056 

-0.4 

Water 

0.492 

986 

-7.0 


Table 3: Calculated Equivalent Permeabilities 


to water properties. Thus, two straight lines are 
observed in the semilog pressure vs. time graph 
[9, 10]. Hence, the reservoir equivalent permeability 
may be estimated using the following equation, 
derived from the single-phase formulation: 

k eq = L1 _ 5lQ: ^ ; ( 40 ) 

hXfnif 

where a p is a unit conversion constant. For a 
consistent set of units, a p — 1, whereas for the set 
of Brazilian field units, a p = 19.03. The straight 
line angular coefficient is denoted by rrif , where the 
subscript / indicates the phase associated to the 
straight line considered. 

Therefore, falloff data were used to generate 
semilog graphs. Pressure data were plotted against 
Horner time. Results are displayed in table 3. 

Equivalent permeability calculated according to 
equation (40) using semilog graphs were very close 
to the input values for all cases, both for water and 
oil properties, expect case Bl. Due to the same 
causes discussed in section 6.3, analytical model 
is subject to noticeable error during early falloff 
times (that means, higher Horner times). Hence, 
the formation of the first semilog straight line may 
be impaired. Case Bl was more sensitive to those 
factors, even though they were present in all cases. 
However, semilog lines associated to oil properties 
allowed the estimation of equivalent permeabilities 
with an error smaller than 1% for all cases. 

It is interesting to notice that, in field tests, 
early falloff data usually present too much noise 
and is subject to storage effect. Thus, water 
semilog straight line may not be detectable, and 


the reservoir equivalent permeability is estimated 
using the straight line relative to oil properties in 
a typical field test. That means the analytical 
data allowed an accurate estimate for the equivalent 
permeability when the most relevant straight line 
(the one associated to oil mobility) is used. 

7. Conclusions 

Based on the formulation for falloff period in 
single-layer reservoirs proposed by Peres et al. 
[12], an expression for multilayer reservoirs was 
developed. The suggested model is analogous to the 
injection period solution, formulated by Barreto et 
al [13]. 

The proposed model was applied on a set 
of cases, with distinct number of layers, layer 
permeabilities and layer skin factors. For all 
cases, comparison between analytical solution and 
numerical simulation showed a close agreement. 

This indicates that, after the well is shut, 
vertical flow between layers through the wellbore is 
negligible, even if layer permeability are different. 
Such strong hypothesis was essential to reach the 
analytical model for falloff period in multilayer 
reservoirs. Results suggest that this fact may not 
be true in reservoirs with significantly different 
individual layer skin factors. 

Total flow-rate along the reservoir were computed 
using different approximations. Falloff pressure 
derivative profiles presented an unexpected result, 
showing that approximating total flow-rate by 
the oil flow-rate yields the most realistic pressure 
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derivative behavior, although this is the roughest 
way to approximate total flow-rate. 

The proposed formulation was also successfully 
used to estimate the reservoir equivalent 
permeability using falloff data and oil properties. 
For all simulated cases, calculated values showed 
an error smaller than 1%. 

Such results verified that the proposed 
formulation was accurate to describe falloff 
pressure behavior in multilayer systems. 

N omenclat ure 

Bf = Phase / formation volume factor 
c t = Total compressibility 
hj = Thickness of layer j 
kj = Absolute permeability in layer j 
k r f = Relative permeability of phase / 

P = Pressure 

P w f = Wellbottom hole pressure 
qj = Total flow-rate in layer j 
qfj = Flow-rate of phase / in layer j 
r = Radius 

rpj = Waterfront radius in layer j 
r w = Wellbore radius 
S w = Water saturation 
t = Time 

A f — Mobility of phase / 

A f — Endpoint mobility of phase / 

<fij = Porosity in layer j 
jifj — Viscosity of phase / 
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